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Abstract 

In present article the self-contained derivation of eigenvalue inverse problem 
results is given by using a discrete approximation of the Schrodinger operator 
on a bounded interval as a finite three-diagonal symmetric Jacobi matrix. This 
derivation is more correct in comparison with previous works which used only 
single-diagonal matrix. It is demonstrated that inverse problem procedure is 
nothing else than well known Gram-Schmidt orthonormalization in Euclidean 
space for special vectors numbered by the space coordinate index. All the results 
of usual inverse problem with continuous coordinate are reobtained by employ- 
ing a limiting procedure, including the Goursat problem - equation in partial 
derivatives for the solutions of the inversion integral equation. 

PACS 02.30.Zz 

1 Introduction 

There is a vast literature on the inverse scattering problem. Suffice it to mention the 
classical monographs [1-5], see also j^j. That theory has multiple applications and 
never stopped developing \J\. In parallel with its renovation, attempts to give a clear 
and obvious treatment were undertaken [8-10]. In these last papers that was done by 
using finite-difference approach which reduces the problem to solving relatively simple 
algebraic equations. Passage to the limit of the continuous variable allows one to obtain 
classical results of the inverse problem. Thus, the finite-difference version represents a 
valuable tool to reproduce all the results of the inversion procedure on a more accessible 
level of understanding. And this is not only of pedagogical interest. The aim of science 
is, among others, to supply a maximally compacted and clear knowledge free from 
superfluous and often obscure details. 

The authors of [8-10] restricted their consideration to the finite-difference matrix 
Hamiltonian with potential coefficients only on the main diagonal of the operator. So, 



there appears a disparity in numbers of interaction parameters and spectral data (see 
discussion below). As a result, we need either to impose some restrictions on spectral 
parameters or to introduce additional non-local potentials as was done in present article. 
That problem was not considered in [8-10], which led to an "erroneous" final result for 
the potential in the continuum limit: there must be additional factor 2 (missed in the 
articles mentioned) in front of derivative of the solution of inverse problem integral 
equation [see Eq. ([37)1] . This oversight was partially compensated in the paper [TT] 
where the authors introduced non-local potentials which are needed for correct final 
results in the continuum limit. However, their method had several deficits, among 
them one can mention a rather "adjustable" character of the procedure suggested and 
none of their algorithms indicates the uniqueness. 

At the same time, the inverse eigenvalue problem for the discrete analog of the 
Sturm-Liouville operator is now well developed. There is sufficiently large number of 
papers on that subject, see ^2] where one can find most of the references. There are 
different variants of the inverse problem in the discrete approximation. All they deal 
with a matrix (finite or infinite) with several diagonals which may differ in number 
and which are recovered by given (spectral) parameters. For the Gel'fand-Levitan 
analog, the most known inversion variant, the central theorem is valid: Given the 
set of eigenvalues of a three-diagonal Jacobi matrix and the first components of the 
associated orthonormal eigenvectors, there exists a unique Jacobi matrix corresponding 
to these data, see e.g. [121 EH and references therein. The usual proof of this theorem 
is performed using orthogonal polynomials [12J. 

Thus, we have the continuous and discrete variants where the inversion procedure is 
well established. However, considerably less was done to link them both. The problem 
is likely that the recovery procedures in discrete and continuum cases outwardly have 
little in common. In the continuum version, the inversion procedure (by Gel'fand- 
Levitan-Marchenko) is built as a transition from a certain known system (free motion, 
as a rule) to the system with known spectral data but with unknown potential to be 
reconstructed. The aim of the paper is to give such a derivation of the inverse eigenvalue 
problem from its discrete variant which would be free from the previous errors. 

In the discrete variant, we have to develop, in a more explicit form, a structure 
similar to that in the continuum case. In doing so, we have to use a general criterion 
which would strongly specify whether our development is correct. The orthogonal poly- 
nomials' method gives such a hint. This is orthonormality relation (in a special spectral 
measure) which is valid for any system. As is well known, the Gram-Schmidt method is 
essential in constructing these polynomials. So the central idea of the present paper is to 
employ that method in order to reconstruct the set of eigenvectors orthonormalized in 
the initial spectral measure when the last is changed in a given way. We shall construct 
some "prototype" of the transformation procedure that realizes the reconstruction of 
the potential of regular Schrodinger operator in the continuum case. 

We begin, in the second section, with several results consisting of some preliminary 
constructions which appear as intermediate steps in the course of the methods dis- 
cussed. First of all, we give the discrete statement of the Sturm-Liouville problem on 
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a bounded interval with zero boundary conditions, which is equivalent to consideration 
of a three-diagonal symmetric Jacobi matrix. In the continuum limit these diagonals 
merge in a single diagonal (local potential). Then we pose the inverse eigenvalue prob- 
lem in terms of eigenvalues and associated spectral weight factors, first components of 
the orthonormalized eigenvectors. Introducing so-called regular solutions admits of the 
explicit presentation of these parameters which serve as a spectral measure (of bounded 
support) entering in the Parseval relation for eigenvectors. That measure allows us to 
represent this equality as ort honor mality condition for the same vectors but from an- 
other standpoint when the energies represent components and the discrete coordinate 
numbers the vectors. Next step is applying Gram-Schmidt technique to obtain the 
new orthonormal set of eigenvectors (in the sense mentioned) corresponding to the new 
measure. We shall see that procedure indeed reproduces prototypes for equations of 
the inverse problem in the continuum limit. The proof is given that such a orthog- 
onalization is the only possible development. Next, we use the new eigenvectors to 
recover the potential coefficients on the three diagonals of the Jacobi matrix (discrete 
Sturm-Liouville operator) by using completeness relation for the new eigenvectors. 

In the third section, we pass to the continuum limit. We demonstrate how all the 
discrete equations-prototypes go over into the classical equations of the inverse Sturm- 
Liouville problem: Gel'fand-Levitan equations, expression for the potential, classical 
Goursat problem, etc. That accomplishes our program. 

2 Discrete version of inverse problem on finite in- 
terval 

It is most easy to demonstrate the essence of the inverse eigenvalue problem for the 
example of finite-difference Schrodinger equation in the discrete variable x n , n G Z with 
the step A: 



V(x n ) and u(x n ) are real, for this problem is reduced to linear algebraic equations. 
The first three terms in this equation represent finite-difference operator of the second 
derivative, i.e. kinetic energy. Note the existence, in the Schrodinger equation, of terms 
u{x n ) corresponding to a "minimally non-local" interaction. We shall soon come back 
to them and their introduction will turn out justified. 

Let us consider the bounded interval [0, ir] with finite number N of points inside: 
xo = 0; xn+i = 7r, so that A = x n+ i — x n = n/{N + 1). Let us add the Eq. (1) by the 
Dirichlet boundary conditions: 




(1) 



*(x ,E) = V(x N+1 ,E) = 0. 



(2) 
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These zero boundary conditions have the well-known physical interpretation that the 
movement of a particle is restricted by the infinitely tall walls at the points xq and xn+i 
(the infinite rectangular potential well). 

The spectrum of the problem (JTJ), (J2J) is a ladder of discrete energy levels {E U }^ =1 
for bound states representing the unit vectors ^ u (x n ) = ^/(x n ,E u ), 

N+l 
ra=0 

The Sturm-Liouville problem (0), (j2J) can be rewritten in a more visible form by 
using the symmetric tridiagonal (Jacobi) matrix (N x N) 

H = f + J, 



T 



( 2/A 2 -1/A 2 
-1/A 2 2/A 2 -1/A 2 
-1/A 2 2/A 2 



\ ■ 

which acts on the vector-column ^ v G 






-1/A 2 



/ V(x l ) u(x x ) 
u(xi) V(x 2 ) u(x 2 ) 
u(x 2 ) V(x 3 ) u(x 3 ) 



-1/A 2 2/A 2 / 



(3) 



u(x N -i) V(x N ) J 



( V(xi,E u ) \ 
y(x 2 ,E t 



This explicit form elucidates why the u(x n ) stands in front of $(x n+ i,£): that is 
because of the upper u-diagonal in J which is one element shorter than the main V- 
diagonal and contains N-l elements. The same coefficients form up the lower diagonal 
(thanks to the matrix symmetry). Let us specially note that homogeneous boundary 
conditions generally different from (J2J) would require some modification of the matrix 
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representation (J3J). In passing to the continuum limit A — > 0, it will be impossible to 
distinguish the u- and V-diagonals, i.e., the resulting interaction will be simply the sum 
of limiting values for u's and Vs. A paragraph later, we shall give the motivation for 
the appearance of the additional diagonals in the interaction matrix J. 

Besides the energy levels, let us introduce additional fundamental spectral parame- 
ters, namely, norming constants or spectral weight factors. By the definition, these are 
the coefficients c v of proportionality between the normalized eigenstates ty u (x n ) and 
the regular solutions <p(x n , E u ) at the eigenvalue energy, ip(0, E) = 0, <p(xi, E) — A [i.e. 
the derivative is equal to 1], ip u (x n ) = ip(x n , E u ): 

^u{x n ) = c u ip u (x n ). (4) 

The continuum analog of the regular solution satisfies ip(0,E) = 0, <p'(x, E)\ X=Q = 1. 
The continuum generalization of the spectral weight factors introduced is obvious. In 
classical inverse Sturm-Liouville problem (continuous coordinate), it is well known that 
the double set of the spectral parameters {E u , c u } uniquely specifies the potential. 

Now we get to the core of our paper. Our task is to pose such a discrete version 
of the inverse Sturm-Liouville problem that, in passing to the continuum limit, as 
direct as possible reproduction of all the result of the continuum version is feasible. 
So, the problem within which it seems logical to work is posed as follows: Given the 
set {E u ,c u } with c u in (jlj), the potential matrix J with coefficients V(x n ) and u(x n ) 
is recovered completely. In principle, a question may arise whether the three-diagonal 
Hamiltonian with the "non-local" u-coefficients is consistent with the uniqueness of the 
potential recovery from the set {E u , c u }. What does force the extra w's? As was already 
mentioned in the introduction, any set {E u , c u } can occur for a unique three-diagonal 
H. Moreover, there is an extension of that result. It is the theorem by Gladwell and 
Willms in which the statement was proved that a symmetric p-band matrix (a 
matrix with 2p + 1 bands, p bands below the diagonal) may be uniquely constructed 
(apart from certain sign ambiguities) from its eigenvalues and the first p components 
of its normalized eigenvectors. Hence, once we know all E u 's and the first eigenvector 
components, by virtue of Eq. these are Ac u in our problem, we can uniquely restore 
1-band, i.e., three-diagonal Hamiltonian (jSJ). That also accounts for the u's. 

Let us give an additional "half-heuristic" explanation of this fact. If we have only 
a local potential V(x n ) with N values at N points, the number N of free parameters 
{V(x n )}n=i equals exactly the number of eigenvalues E u . To the point, the correspond- 
ing inverse problem has not a complete solution. If we introduce both spectral param- 
eters, E v and c u , while the single diagonal J persists, we shall really face a problem 
of over-determination of the set {E u , c u } that contains 2 N-l free parameters. In fact, 
there are N levels E u and N-l parameters c u by virtue of the relation Y^v=i c u = 1/A 3 
that follows from © for n — m — 1, while V(x n ) does only N ones. It is introduction of 
N-l coefficients u(x n ),n = 1, ...N — 1 into Eq. ((TJ) [or additional diagonals in 0] that 
ensures the equality of numbers of spectral data and interaction parameters. In the case 
of continuous coordinate, overfilling the set of spectral parameters reveals itself only 
in many-dimensional D > 2 problems, so we shall manage to restore one- dimensional 
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local potential by the complete spectral set [see (|57j)]. 

The functions (p v (x n ) can be considered as vectors in special Hilbert (Euclidean, 
to be precise) space, in which the coordinate x n numbers the eigenvectors and energy 
index v is only used to denote the z/th vector component. The inner product in that 
space is determined by the measure given by the spectral weight factors c v . In fact, the 
Parseval's completeness relation 

N 

J2 Vu(x m )V v (x n ) = 5 mn /A (5) 
can be rewritten using Eq.(jU) as 

N 

y^dt(p v (x n )<p u (x n ) = 6 mn /A. (6) 

v=\ 

Let us consider this expression as an orthogonality relation for the vectors (p u (x m ) and 
^(xn) (in the limit A —>■ 0, the "numbers" x m>n of the vectors become continuous 
variable x). Here, the inner product is given by not simply a sum over energy index v 
but a sum with a weight (measure) c^. 

Different potentials correspond to different weight factors determining the metrics 
of our "energy space" but the relation (0) holds true for any potential. In the classical 
variant, the inverse problem can be treated as a transition to the sought potential 

o 

V (x n ) V(x n ) from a certain "initial" (in what follows we shall use the symbol "o" 

o 

to denote everything related to the initial system) potential V (x n ), for which all the 

o ° o __ 

solutions <f u (x n ) and the whole spectral set {E u ,c u } are known, and the relation (p) 
is valid: 

N 

° c ^» = <w a - (7) 

u=l 

All this gives us a hint for deriving new solutions corresponding to the given spectral 
set {E u ,c u }^ =l . Although we do not know yet the sought potential matrix J, we 
beforehand know that the regular solutions (p u (x n ) to these potentials must satisfy 
the orthogonality relation (JOJ) with the new c„. We shall catch at this fact and use 
the orthogonality relation as a central criterion in finding new eigenvectors (solutions 
^(in)). Changing the metrics of Euclidean space in replacing c v — > c u results in that 

o 

the "old" unit vectors <fi u (x n ) are no longer orthogonal. So the idea is as follows: Once 
new unit vectors must satisfy Eq.(fS|). we could obtain them, e.g., orthogonalizing the 

O 

¥v (x n ) with the new weight (? v by the Gram-Schmidt scheme. In other words, the new 
vectors obtained by that way and satisfying (JHJ) with the weight multipliers <? v will be 
the solutions to the new potentials V^'s and w's. Indubitably, this makes sense only 
when the procedure really gives the desired vectors, i.e., it is unique (see the proof 
further on). 
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For simplicity, we shall at first think E u = E v . Let us recall this standard orthog- 
onalization procedure for the example of two initially non-orthogonal (in sense of new 

o o 

weight function) vectors (i.e. when N=2) [x\),^P v {x 2 ),v = 1,2. As a first unit 

o 

vector (p v (xi) of the new system, we take the unchanged unit vector <p u (xi), and the 
second unit vector is constructed from the second unaltered one, only we have to sub- 

o 

tract everything superfluous (parallel to <P V { x l))i f° r the orthogonality with the new 
measure: 

00 o 

Vv(x x ) =V V {xx)\ ip v {x 2 ) =¥„ (x 2 ) + AK(x 2 , xi) V u (xi). 

The coefficient K(x 2 ,Xi) is derived from the condition of orthogonality of the new 
vectors with the new weight c u : 

o 

(p v (x 2 ) -L ifuiXx) =fp v (xi). 

We have 

2 2 o2 

°l ^ ( x2 ) + A Yl <% K ( x 2> x i) m = 

u=l u=l 

2 2 2 o2 

K(x 2 , Xl ) +^4<P V (xi) <P V (x 2 ) + AK(x 2 , Xl ) ^(c 2 - cj <f v (xx) = 0, (8) 

o2 

where we add and subtract the term c v from the multiplier c v and, furthermore, use 
Eq. (|7jl. We can rewrite the last equality in the form as follows (extremely simplified 
two-dimensional "prototype" of the inverse problem equation): 

K(x 2 , xi) + Q(x 2 , Xl ) + AK(x 2 , x 1 )Q(x 1 ,x 1 ) = 0, (9) 



where 



E2 ° ° 
£p (•I'm) ^Pv 



N=2 

9 / n 
(•'•„)- 

i/=l 

5^ (x m ) ^ (x„); m, n = 1, 2. (10) 

/Lt=l 



In general case of N-dimensional Euclidean space we shall follow the same scheme. 

o 

In doing so, it is possible to take into account the case when the levels change: E u ^ E v . 



We must orthogonalize N vectors by the measure c 2 u : <-P (xm^E^), (m = 1, 2, ...N, u = 
1, 2, ...N). Consequently, we have for the new solutions 



m— 1 



tp(x m ,E v ) =£> [x m ,E v ) + ^AK(x mi x n ) [x n ,E v ), (11) 



n=l 
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where the coefficients K [the kernel of the transformation operator (jllj) ] follow from 
the conditions of the orthogonality of new vectors (by measure cl) ip(x m , E v ) (m = 
1,2, ...iV): 

ra>ni E v ) J_ ip(x n ,E u ), 

which lead to the system of algebraic equations for K - discrete analog of central 
equations of the inverse problem: 

m— 1 

K(x m ,x n ) +Q(x m ,x n ) + ^2 AK(x m ,x p )Q(x p ,x n ) = 0, m > n, (12) 

P =i 

where Q(x n ,x m ) is determined as in (JTUJ), only the values m and n are no longer 
restricted by 1 and 2, and the indices \i and v number solutions at initial and shifted 
energy levels, respectively 

N 



E2 



N 



° c l$ {Xm,E^) <P {x n ,E^). 



(13) 



Let us note that the form (fTTj) ensures the desired boundary condition for the regular 
solution: <p(xo, E v ) = 0; <p(xi, E v ) = A. The system Eq. (|T2j) of the recurrence compu- 
tation of the K's provides them uniquely. We can formally introduce the diagonal terms 
K(x n ,x n ) (bearing no relation to ip(x,E v )) such as K(x n+ i,x n ) — K(x ni x n ) ~ 0(A) 
which will be useful in what follows. 



When the first m + 1 unit vectors <£> (xj, E v 



1, .., m + 1 are orthogonalized, this 



corresponds to an intermediate submatrix-block transformation of the initial Jacobi-like 
operator (JHJ) so that 



J 



J r, 







Jtv- 



•J [> 



( V[x x ) u{xi) 
u{x\) V(x 2 ) u(x 2 ) 
u(x 2 ) V(x 3 ) u(x 3 ) 



\ 



o 

Jn- 



\ ■ .0 u(x m _x) V{x m ) u (x m ) ) 

f U (x m ) V (Xm+l) U (x m+ i) N 

u (x m+ i) V {x m+2 ) u (x m+2 ) 



u (x N -i) V (x N ) J 



(14) 
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where the symbols in the top-right and bottom-left corners of the J-matrix denote 
zero (m x N — m — I) and (N — m x m — 1) matrices, respectively. The submatrix 

O 

J m is formed up by the perturbed coefficients while the Jn-tu is not affected yet by 
the transformation associated with the reorthogonalization. Note that the last row of 
the submatrix J m contains only two transformed elements, u(x m -i) and V(x m ), apart 
from the u (x m ). This is because the (m x m) quadratic submatrix was transformed 
only which contains, in its last row, two non-zero elements mentioned. The coefficients 
of that intermediate transformation block may be found from formulas (J2Uj) and (|21|) 
where one should substitute u (x m ) for u(x m ). See also the formulas (|2Uj) and (J2JJ) and 
the subsequent discussion. 

Now we are ready to give the proof that the above procedure is unique. Let us carry 
out it by induction. Suppose that the desired transformed vectors are uniquely given 
by the equation (jllj) for m < N, N < N for a certain N being the integer. For N = 1, 
this is verified trivially. Let us show that the formula (fTTj) holds true at the point x^+i- 
Indeed, the (p(x^ +l , E u ), being orthogonal to all <p(x m ,E v ),m < N, can be sought, in 

o 

principle, as a combination of the initial V 9 (x n , E v ) for all n. The coefficients of such 
a hypothetical combination (there are N pieces of them in all) have to be determined 
from the condition of orthogonality of (p(x^ +1 , E u ) to both <p(x m , E v ),m = 1, ...N and 
certain N — N unknown vectors from the new orthogonal set For ascertaining 

what we shall do now, we involve the Schrodinger equation (JJJ, which is a recurrence 
procedure of the step-by-step computation of the <p(x n ,E v ). Specifically, we mean the 
block operator in ([Ti|). last row: 

<p{xfr +1 , E v ) - 2y?(xjy, E v ) + y?(sjy_i, E v ) 

^2 *~ v ^nip^n, A/ j 

+ u {xn)(p(xn +1 , E v ) + u(xn-i)^(xn-i, E v ) = E v if{xfj, E v ). 

It is seen that u (xft)(p(xft +1 , E v ) is a linear combination of u(x^_i)(p(x^_i, E v ), 
V(xft)(p(xx, E v ) and E u ip(x^, E v ). Consequently, by the assumptions of the validity of 

— o ° 

(fTTj) for m < N,u (xft)(p(xx + i, E v ) can only be represented through V (x K [x # , x m ) V 9 

— o o 

(x m ,E u ), m = 1, ...N - 1; V(xfj) V 9 (xft,E u ); u(xff_ 1 )K(xff_ 1 ,x m ) ¥ (x m ,E„), m = 

— O O — 

l,..JV-2; w(xjv-i) V 9 {xn-i,E v ); K(xf I _ ll x m )E v V 9 (x m ,E u ), m = 1,...N-1; 

o 

and E v V 9 (xff,E„). For the last term, we find from the non-perturbed Schrodinger 

o o ° ° ° 

equation that E v V 9 {xfj,E v ) is expressed through u (x$) V 9 (xn+i,E u ), V V 9 

° ° ° 

(xjv, E v ), u (xjv-i) V 9 (£jv-i) E u ) [and similarly for other E v <P [x m , E v )\ and, finally, the 

O — 

<f(xft + i, E v ) must be sought as a linear combination of V 9 (x m ,E u ), m = 1,...N + 1. 
In other words, <p(xft +1 , E v ) E spanjy 9 (x m , E u )}^t\, besides that <p(xft + i, E v ) _L 

o . 

spanjv 9 (x m , E u )}^ =1 = span{y9(x m , ^)}„ =1 . It is well known that the Gram-Schmidt 
orthogonalization enables a unique solution satisfying these two conditions. Note that 

— o o ° 

(_/V+l)th term u (xjy)<p(xft +1 , E v ) is a combination of the summand m (xn) V 9 (xx +1 ,E u ) 
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(with the same coefficient) and other terms with m < N. In other words, in the de- 

o 

composition of V'wv+Ij E u ), we have the term V 9 (x^ + i,E u ) with the unit coefficient. 
Consequently, <f(x^ +1 , E u ) is represented in the form (|TT|) again (QED). 

It should be noted that the formula (fTT|) is also valid for solutions at energies E lying 
between the levels E v where the regular solutions, being the Cauchy problem solutions, 
are defined (though non-physical). Do not confuse "running" energy values at which 
the solutions ip(x m ,E) are defined with the energies occurring in the inverse problem 
equations (fT3|) . Indeed, let us decompose (f(x m , E) into the complete set of the solutions 
<p(x m , E u ) (in sense of the usual inner product ^^ =1 Ac u c^ip(x m , E u )<p(x m , E^)): 

N N 

(p(x m ,E 1/ )(p(x m ,E). (15) 

u=l m=l 

Since the ip(x m , E u ) J s are expressed, in accordance to (fTTj). through the unperturbed 

o 

f (x m , E v ) then we shall expand them, too, in a complete set of the old solutions 

o o 
[p^mi Eft)- 

N N 

<P (x m , E v ) = C(E V , E^) V (x m , E^); C(E V , £ M ) = ^ A 

C'fCp y^rn: E u ) V {x m ^Efi)- 

fi=l m=l 

(16) 

Combining |TT| , (|15jl and (|T6|) , we get the following expression for the new solutions 
at arbitrary E: 

m—l 

ip(x m , E) = A ( E , E u, E^) i {x m , Ef,) + J2Y1 AK ( X ™, x n ) 

fi,u ft,v n=l 

xA(E,E v ,Ev)$>(x n ,Ev); A(E, E u , Eft) = £(E, E U )((E U , Eft)- (17) 

In the limit when the new and old spectral parameters coincide, K vanishes and, hence, 
<f(x m , E) turn into unperturbed solution 



¥ (x m , E) = J2 A(E, E v , E^) £ (x m , Eft) 



fl,V 



Substituting this expression in (fTTj) . we get: 



771—1 



<p(x m , E) =tp (x m , E) + J2 AK (xm, x n ) V (x n , E). (18) 



71=1 



Let us stress here that K is independent of energy E. The formulas ()18j) . f!12l) and 
(HHJ) give the expression for K in the form of sum of products of the old solutions and 
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transformed ones: 



N 



N 



K(Xr 



, Xr 



(19) 

It remains to obtain equations for the transformed potentials V and u. We already 
know the solutions of Eq. with the unknown potentials V(x n ) and u(x n ) [see formu- 
las (|TT|) and (fI2|)]. i.e. eigenvectors of the new Hamiltonian plus associated eigenvalues 
E v . As was mentioned, by virtue of the theorem by Gladwell and Willms |13| . that 
is enough for the three-diagonal Hamiltonian matrix (with off-diagonal elements) to 
be uniquely recovered. These authors used the block Lanczos algorithm. However, we 
shall apply an outwardly different method pursuing the aim of reproducing final for- 
mulas in the continuum limit. Let us multiply both parts of the Schrodinger equation 

o o o o o o 

(HJ) for the solutions ip(x m , E^) [Eq. (JTHJ)] and V 9 (x n , E^) by V 9 (x n , E^) and (p(x m , i? M ), 

respectively, sum over /i with weight c a and subtract from each other the resulting 
expressions. At fixed m, we perform this procedure for n — to, to — 1, .... In calculating 
sums (over /i) one should take into account the relation ((7j). As a result, we get to the 
following equations for V and u: 

o 
{V(x m )- V {x n )}K(x m , X n ) + u(x m )K(x m+ i, X n )- U (x n ) 

xK(x m ,X n+ i) + u(Xm-l)K(Xm-l, X n )— U (x n -l)K(x 

mi Xn—l) 

K-(Xm+li Xn) 2i^(x m , X n ) ~\~ K(x m — 1, X n ) 



A 2 

K-{Xrm X n +l) *2K (Xfyii Xn) ~\~ K(x m i X n — l) 



A 2 



n < m — 2: 



(20) 



and for n = m, m — 1 



)-u(x m _ 1 ) _ K(x m+ i,x m - 1 )-K(x rn ,x m - 2 ) 

A A 2 

o 

-V(x m )K(x m , X m -l)+ V (x m -l)K(x 

mi Xfn—% ) 

-u(x m )K(x m+ i,x m -i)+ u (x m ~2)K(x 

mi X m —2) i n> 



V(x m )-V{Xm) 

A 



K (^-m+1 t^m) K\X.fn i%m — l ) 

A 2 



~ u (x m )K(x m+ i, x m ) + K(x m , x m -i) u (x m _i), n 



m 



TO, 



(21) 



where the terms K(x m , x n ) for which to, n > N or to, n < 1 are omitted. But for 
n = to + 1 we obtain that w(x m ) —u (x m ). There is nothing strange in it because 

o o o o 

summation is carried out for the term <p (x m +x, E^) that is orthogonal to all <fi (x n , Ep), 
n < to + 1. But the kernel K containing all the information about new solutions just 
stands at these summands, see formula ()18|) for n = to + 1. In other words, for the case 
n = to + 1 the summation expunges everything that bears a relation to the new system 
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under construction. Indeed, from the Schrodinger equation for & (x m+ \, E^) multiplied 

o Q 2 

by <f(x m , Efi) and summed over fi with the weight we have 

N 

E Q 2 o o o o 

C M £ M if{x m ,E^) V (X m +1,E^ 

How should one treat this? One variant is that u(x m ) =u (x m ), which corresponds to 
restoration ofamxm submatrix for which the element u (x m ) is an outer one, see (|14j) . 
This can also serve as a proof of the ()14|). Second interpretation is that u(x rn ) =£u (x rn ) 
(but not for m = N), nevertheless (i.e, the above procedure does not work). This takes 
the place when we have the whole J-matrix transformed but the equations (|21|) by 
themselves did not allow the computation of u(x m ). How then to uniquely restore the 
J will be discussed a bit later, but now we must ascertain for ourselves how it is possible 
that the same solution (p(x m , E) may satisfy the Schrodinger equation (JTJ) with different 
potential coefficients u (x m ) and u{x m ) at x m+ \. The matter is that we deal with (finite- 
difference) non-local potential and this ambiguity is just characteristic of it. Indeed, let 
ip(x m ,E) satisfy the Schrodinger equations with both {ui(x m -i),Vi(x m ),ui(x m )} and 
{u2{x m -i), V2(x m ), U2(x m )}. Subtracting these equations from each other we have 

[Vi(x m ) - V 2 (x m )](p(x m , E) + \ui(x m ) - u 2 (x m )]ip(x m+1 , E) 
+[«i(ar m _i) - u 2 (x m ^i)](p(x m _ 1 , E) = 0. 

If we have only one (local) potential coefficient then it would be the same. But now, 
for several non-local potential coefficients coupling neighbour x-points, that equation 
clearly demonstrates that Vi(x m ) — V2(x m ) and others may all be non-zero. Summing 
up this discussion, we have elucidated that the procedure used for derivation of (|2Uj) and 
(|21|) cannot distinguish all the variants of w(x m )-coefficient determination proceeding 
from a general incapability of giving a unique non-local interaction associated with a 
certain solution of Schrodinger equation. 

However, for the whole vector (p(x m ,E v ), i.e. the solution defined at all the points 
x m , m = 1,...,N we are able to uniquely derived the quadratic potential matrix J 
whose eigenvectors are <p(x m ,E v ). Taking m = N we first find u(a;jv-i), V(xjy) and 
u(xn). Of course, this requires the knowledge of u(xn). But we have no more equations 
for determining the potential coefficient u(xn)- However, we see that the u(xn) is a 
continuation of the last iVth row of the matrix J. This resembles the case with the 
unfinished restoration of the J, see f|14j) . i.e. the potential perturbation (in form of a 
quadratic matrix) never reached the u(xn). That is, the u(xn) is independent of the 

° o 

transformation generated by i^-coefficients. But then, taking {E u ,c u } = {E u , c u}, we 
see that K — and u(xn) exactly corresponds to the reference potential. Thus, we 
have u{xn) =u {x^). Next, at the point xn we have, instead of (|21|). the following 



u (x m ) ^ J_ 
A A 3 ' 
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system of equations 



f at- o2 o o o 

n(xAf_i)-n(xAf-i) _ Y^L i =i c n ( P( x N+i,EL i )¥(x N „ 1 ,E fi )-K(x N ,x N - 2 ) 




o 



-\/(x 7V ) J R'(x A r,a; A r_ 1 )+ F (x N - 1 )K(x N , X N _i) 




k -«(a;jv) Z^=i C M ^(^at+i, V 5 (xjv, E^) + K(x N , xn-!) u (xjv-i), 



(22) 



o 



where <p(xn+i, Efj) is found from the Schrodinger equation ((H): 



o A 2 ° 
<p{x N+ i,Eft) = t~2 7 rNzjv-iMzjv-ijE'/J 




o o o 

+V(x n )ip(x N ,E^)- E^ <p(x N ,Ep)] 



(23) 



o 



<f(x n ,Efj,) given by (fTHj) and u{xn) =u (x^). From (j2~2~|) and (|2^j) we obtain V(xn) 
and m(xat_i). We then substitute the value u(xn~i) (by virtue of the symmetry of 
potential matrix) into equations (|2*U|) for m = N — l;n = m — 2, m — 3 from which 
we find, in turn, V(xn-i) and u{xn-i)- Afterward, we substitute this last coefficient 
into equation (|2*UJl for m = N — 2;n = m — 2,m — 3 and get V(x^^ 2 ) and u(xn- 3 ) 
and so on. Thus these equations allow the computation of V and u via the solutions of 
the inverse problem equation (|T2*|) - the coefficients K(x m ,x n ),m > n (plus additional 
requirement concerning u{x^)). For any finite N, these linear equations are uniquely 
solved. But, with the N large, the numerical instability increases that leads to the well 
known problem of the ill-posed inversion procedure in the continuum limit. However, 
let us drop the discussion on that problem here, especially as we only want to reproduce 
the expression for the continuous potential. So we shall keep on dealing with equations 
(|2()j) and (|21jl and next show that passing to the continuum limit in these equations will 
lead us to the classical results of the Sturm-Liouville inverse problem. 

3 Continuum limit 

Let us now pass to the limit of the continuous variable x, i.e. to the limit A — > 
(N -> oo) so that AN = nN/(N + 1) -> n in the formulas (fUJ), JHJ), (J2DJ) and (j2Tj) . 
Let us recall the standard rules of the transitions from the finite-difference operators to 
their continuum counterparts: 




(24) 
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f{Xn) ~ f{Xn-l) . df ^ 



A dx 
f(x n+1 ) - 2f(x n ) + f{x n ^) d 2 f 



(26) 



A 2 dx 2 ' 

Now let us look at the Parseval's relation that takes, in the continuum limit, its 
usual form for the infinite-dimensional (Hilbert) space 



00 n 
a 2 o 



E„) V (y, E„) =6(x-y), (27) 
n=i 

and the same is for the new regular solutions <p(x, E). 

In the continuum case, we have ip(0,E) = 0, ip'(0,E) = 1. Spectral weight factors 
are in that case, too, the coefficients of proportionality between normalized eigenfunc- 
tions and regular solutions. That is why they are also referred to as norming constants 
since the multiplication by c u = 1/ J* (p 2 (x, E u )dx turns regular solution (at E = E v ) 
into the normalized one, 

*'(x,^)| x =o = c v . 

The expression for the transformed regular solutions has now the following form 
(using 



<p(x, E v ) =V (2, E v ) + / K(x, y) V (y, E u )dy } (28) 

and similarly 

<p(x, E) (x, E) + / K{x, y) <P (y, E)dy, (29) 



where x G [0, 7r]. These formulas have just demonstrated that the passage to the limit 
A — > does exist. For the kernel K of the operator (|28j) which transforms the solutions 
to the initial potential into the solutions to the new one (generalized shift operator), 
we have the continuum analog of Eq. (fT2*|) - the inverse problem equation proper: 

K(x,y) + Q(x,y)+ K{x, z)Q(z,y)dz = 0, (30) 
Jo 

where the kernel Q is constructed from the unperturbed functions with the old and new 
spectral parameters [as in Eq. (|TU|)]: 

Q(x, y) = J2 c l$ (*> E ») * E ») ~ ° c l^ * (y> • ( 31 ) 

For the continuous coordinate, the expression (|T9*j) for K has a similar form: 

K(x, y) = -J2 c ^(x, E v ) V {y, E u ) + J2° c l <P( X > E») V {y, Ej. (32) 
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As we have just carried out the passage to the continuum limit in solutions (fTT|). it 
is clear that such a limit exists for the potential, too. In fact, the expressions for the 
potential coefficients (|2U|l and (|7T|) were secondary with respect to (fTTj). i.e. we always 
'extract' them from Schrodinger equation using the information about its solutions 
(see the section above). As was shown, this procedure is substantially based upon 
the completeness relation which stands good for any A, including the continuum case. 
Moreover, we might use the continuum solution an d the Parseval's relation ()27|) 
for the continuum potential to be derived. However, we choose the way of continuum 
passage in (J20|) and (|2Tf . Another point is that the continuum potential is local. Indeed, 
the u's and V in each row of the discrete Sturm-Liouville operator are specified at the 
very neighbour points x n and x n+ \ merging if we pass to the continuum limit, which 
entails, in turn, superimposing the potential coefficients at one point: V + 2u. The 
distinct feature of the local potential is that the limiting equation must determine it 
uniquely in contrast to (J2Uj) and (j7T|) which, by themselves, could not uniquely specify 
u's and V's by reason of a 'non-local' character of the interaction in the discrete case (we 
remember that there was required an additional knowledge of u(xn) at the boundary 
of the interval for the uniqueness). But for the Schrodinger equation with a local 
potential, it is well known that the potential always occur for the unique solution (with 
given boundary conditions) and vice versa. Thus, we can beforehand anticipate an 
expression for a unique specification of the local potential in the continuum case. 

We shall now prove that, as A — > 0, the equations (J20j) and PT|) go over, respectively, 
into 

{V d (x)- V d (y) + 2[u d (x)- u d (y)]}K(x,y) 

d 2 d 2 

;K(x,y)-—K(x,y), (33) 



dx 2 dy 1 



and 



V d (x)- Vd (x) + u d {x)- u d (x) = 2£K(x, x) 
{V d (x)~ Vd (z) + u d {x)- u d (x)}K(x, x) ( 34 ) 
= ~q^K{x, y) \y= x — Qy-2K(x,y)\ y=x , 

where V d (x) = lim A -»o V(x m ) and, analogously, u d (x) = lim A -»o u(x m ). The tilde sign 
stands for the potentials obtained in passing to the limit of continuous coordinate in 
the solutions of Eqs. (J5TJ). 

In developing these equalities, it is useful to employ the diagonal terms K(x n , x n ) 
such as K(x n+ i, x n ) — K(x n ,x n ) ~ 0(A). First of all, let us consider the term 
(K(x m+ i, x m -i) — K(x m ,x m -2))/ A 2 in (J2H). We add to and subtract from the ex- 
pression in the numerator the term K(x m ^i, z m _i) — K(x m , x m ). Then 

■K\p£m+li Xm—l) ■^■y^mi -Em— 2) ^(^m+l) ^m-l) "I - ^(^m-lj ^m-l) 



A 2 A 2 

K-(x m , Xm) -\- K {x m —\^ x m _i) K(x m , x m ) -|- K{x m ^ X m —2) 

X 2 
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Next, let us again add to and subtract from the new expression in the numerator the 
term 2K{x m ,x m - 1 ): 



A 2 

■K(Xmi Xjnj 2,K (x m , X m _i) -|- K{x m ^ X m — t i) 



A 2 

a 2 



The first two lines in this expression are the second derivatives with respect to the first 
and second argument of K(x,y) (see ()26p ). Hence, in the continuum limit they become 

d 2 d 2 



K(x,y)\ y=x - —K(x,y)\ l 



The third fraction diverges as A — > 0: A~ 1 dK(x,x)/dx. As a result we have 

K(x m +i, x m -i) ~ K(x m , Xm-2) d 2 K( xi 92 K( xi 

^2 dx 2 y>\v= x dy 2 y>\y= x 



+A- 1 ^-K(x,x). (35) 
dx 



Likewise, it is not difficult to obtain that 

■K{Xm+li %m) K{Xmi X m —l) 

A 2 



A 2 



dr x -^K(x,x\ A->0. (36) 
dx 



In equation (|2()jl we see the finite-difference second derivative in an explicit form. So 
in continuum case this equation becomes If we introduce V(x) = V^(x) + 2ud(x) 

o 

the term in front of K(x,y) is simply the difference V(a;) — (y). It is obvious that 
we introduced a local limiting potential which results from the limiting merging of 
V-diagonal and nearby u-diagonals. 

Now let us multiply both sides of the equations of (|21j) for n = m — 1 and n = m 
by A. We sum the resulting equations and pass to the continuum limit. Then, by 
virtue of (|33j) and (j3T?j) . we get to the first equation in (jMj) valid to within 0(A) (the 
multiplication by A has removed the divergence associated with A -1 ). 

The last equation in (J34)) is not obvious. Indeed, one would think that the term 
u(x m ) must first be derived from the recurrence procedure (|2U|) and (J2TJ) and only 
afterwards can the passage to the limit A — > be carried out - the procedure of a 
prodigious complexity. However, we find a way out: we simply take the sum of non- 
diverging terms (taking into account the expression (|35Jl ) in the right-hand side of the 
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equation (|21|). n = m — 1 to be zero in the continuum limit, i.e., we get the last equation 
in (J211). This by no means contradicts the uniqueness of the sought limiting potential. 
First, this provides the limiting (for x = y) equation for K(x,y) which must exist, 

° o 

obviously. Second, by continuity, the factor Vd{x) — Vd (y) + 2[ud{x) — u d (y)\ in front 
of K(x,y) must coincide with Vd{x) — Vd (x) + Ud{x)— u d (x) when x = y. In other 
words, that means that V(x) = Vd{x) + 2ud{x) = Vd(x) + Ud{x)+ u d (x), i.e. the 
solutions of (f2"T]) go over, in the limit A — > 0, into the same local potential V(x), which 
was beforehand clear. Hence, with the new definition for V(x), we have from (}3*3*j) and 
(El: 

V(x) =V (x) +2—K(x,x). (37) 
ax 

This is the known result of recovering potential in continuum case, which only now 
became reproducible from a discrete mathematics. 
The Eqs. ()33|) and ()34|) can now be rewritten as 



{V{x)-V(y)}K(x,y) 

Q x 2 K {X ', ?/) 

) J 

' dx 



£,K(x,y)-^K(x,y) (38) 
V(x)-V(x)=2±K(x,x) 



This system (added by K(0, 0) = 0) represents the classical Goursat problem [for de- 
termining K(x, y)} and its solvability follows from well known theorems. 

The orthogonalization can also be started from the last vector <p(7r) with 'number' 
x = tt at the right boundary of the interval [0,7r]. Then, instead of solutions p(x), 
the solutions f(x) will be used such that f(n) = 0, /'(7r) = 1. The corresponding 
inverse problem equations, that can be associated with the orthogonalization 'from the 
right to the left', have analogous form as Eqs. (j^Hj) . (|3Uj). (|3*Tj) and(|3*7j). only with other 
integration limits and different sign in front of the derivative in the expression for V(x): 

f(x, E) =°f (x, E) + f K(x, y) f (y, E)dy; (39) 

J X 



K(x,y) + Q(x,y)+ K{x, z)Q(z,y)dz = 0; (40) 



Q(x, y) = Y,ll°f (x, E v ) f (y, E v ) - V (x, E„) f (y, E„)\ (41) 



V(x) =V (x) - 2—K(x, x). (42) 
dx 
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Here, the symbol 7„ stands for the spectral weight factor which is analog of c v . The 
only discrepancy is that the 7„ characterizes the behaviour of eigenfunction at the right 
boundary: 

V{x,E v )= lv f{x,E„), 7, = ^(x,E u )\ x=w . 

At last, let us mention about the eigenvalue inverse problem for Schrodinger equation 
added by boundary conditions of arbitrary kind: 

tf'(O) - (0) = 0, *'(tt) + G^(tt) = 0. (43) 

Here we also have analogous inversion equations, and as spectral weight factors there 
appear the values of corresponding eigenfunctions at the interval edges: c u = ^(0,E U ) 
or 7v = *(ir,E v ). 

4 Conclusions 

In the present paper we carried out the derivation of main formulas of the inverse 
eigenvalue problem on the base of its discrete approximation. Several statements of 
that problem are developed by now, we selected such a statement in which it is possible 
to reproduce in a maximally straightforward way the future structure of the limiting 
inversion procedure: the transition from a known system to the system with given 
spectral data (eigenvalues plus norming constants) but with unknown potential to be 
recovered. The off-diagonal elements are introduced into the matrix Sturm-Liouville 
operator (three-diagonal matrix), which is consistent (in contrast to previous works) 
with the problem statement involving this double set of spectral parameters. At last, 
in comparison with usual derivation of the continuum inversion equations, our devel- 
opment seems to be none the more complicated. At the same time, the reader acquires 
the ability to track in more detail additional aspects of the formalism, in particular 
to look upon the operator transformation realizing the recovering procedure as the 
orthonormalization of the operator eigenvectors. 

References 

[1] Agranovich Z S and Marchenko V A 1963 The Inverse Problem of Scattering 
Theory (London: Gordon and Breach) 

[2] Levitan B M and Sargsjan I S 1975 Introduction to Spectral Theory (Providence, 
RY: Amer. Math. Soc.) 

[3] Marchenko V A 1986 Sturm- Liuoville Operators and Applications (Basel: 
Birkhauser) 

[4] Chadan K and Sabatier P 1989 Inverse Problems in Quantum Scattering Theory 
2nd edn (Heidelberg: Springer) 



18 



[5] Newton R G 1982 Scattering Theory of Waves and Particles 2nd edn (NY: 
Springer) . 

[6] Freiling G and Yurko V 2001 Inverse Sturm- Liouville Problems and their Appli- 
cations (NY: NOVA Science Publ. Inc.) 

[7] Sabatier P S 2000 Past and future of inverse problems J. Math. Phys. 41 4082-124 

[8] Case K M and Kac M 1973 A discrete version of the inverse scattering problem 
J. Math. Phys. 14 594-603 

[9] Case K M 1973 On discrete inverse scattering problems. II J. Math. Phys. 14 
916-20 

[10] Case K M 1974 Orthogonal polynomials from the viewpoint of scattering theory 
J. Math. Phys 15 2166-74. 

[11] Berryman J G and Green R R 1978 Discrete inverse scattering theory and the 
continuum limit Phys. Lett 65A 13-15. 

[12] Gesztesy F and Simon B 1997 M-functions and inverse spectral analysis for finite 
and semi-infinite Jacobi matrices J. d 'Anal. Math. 73 267-297. 

[13] Gladwell G M L and Willms N B 1989 A discrete Gel'fand- Levitan method for 
band-matrix inverse eigenvalue problems Inverse Problems 5 165-179. 



19 



